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Abstract 

The nonextensive statistics based on Tsallis entropy have been so 
far used for the systems composed of subsystems having same q. The 
applicability of this statistics to the systems with different q's is still 
a matter of investigation. The actual difficulty is that the class of 
systems to which the theory has been applied is limited by the usual 
nonadditivity rule of Tsallis entropy which, in reality, has been estab- 
lished for the systems having same q value. In this paper, we propose 
a more general nonadditivity rule for Tsallis entropy. This rule, as the 
usual one for same g-systems, can be proved to lead uniquely to Tsallis 
entropy in the context of systems containing different ^-subsystems. 
A zeroth law of thermodynamics is established between different q- 
systems on the basis of this new nonadditivity. 

PACS : 02.50.-r; 05.20.-y; 05.70.-a 

1 Introduction 



s q = ^= lF \ (1) 



The starting point of the nonextensive statistical mechanics (NSM) of Tsallis[l] 
is the entropy given by 

l-q 

where the physical states are labelled byi = l,2,...,iw,gisa parameter char- 
acterizing the nonextensivity of the theory, and pi is the probability for the 
system to be found at state %. The application of the principle of maximum 



1 



entropy under appropriate constraints can lead to power law distributions 
characterized by the so called g-exponential functional [ ]. For example, for 
canonical ensemble under the constraints associated with probability normal- 
ization and average energy, we can have [2, 3]: 

Pi = -^-[i-(i-?)W 1/1 - fl N>o (2) 

where Z q is the partition function, Ei is the energy of a microstate i, and (3 
is an intensive variable analogous to the inverse temperature in Boltzmann- 
Gibbs statistics (BGS). j3 is the central quantity discussed in this work. 

For a composite system (A+B) whose joint probability pi(A + B) is given 
by the product of the probabilities of its subsystems having the same q, i.e., 
Pij(A + B) = Pi(A)pj(B), we have, from Eq.(l) 

S q (A + B) = S q (A) + S q (B) + (1 - q)S q (A)S q (B) (3) 

and, from Eq.(2) 

Eij(A + B) = Ei(A) + Ej{B) + (q- l^E^E^B) (4) 

Eqs.(3) and (4) characterize the nonadditivity of the same g-systems having 
Tsallis entropy. In other words, they prescribe the class of nonextensive 
systems to which NSM may be applied. They should be considered as the 
analogs (or generalizations) of the entropy and energy additivity in BGS and 
play the role of starting hypothesis specifying the limit of the validity range 
of NSM. It has been proved that [4, 5] Eq.(3) uniquely determines Tsallis 
entropy and that [6, 7] Eqs.(3) and (4) are a group of necessary conditions for 
the existence of thermodynamic stationarity between nonextensive systems. 
When q — 1, S q — Si — — J2T=iPi^ n Pi^ an d Eqs.(3) and (4) become the usual 
additivity assumption of entropy and energy for BGS. 

The range of validity of this generalized theory is actually under study. In 
view of its application to complex systems like turbulence, economic systems, 
fractal and chaotic systems, a point of view [8, 9, 10] is more and more 
accepted that NSM can be used to describe some nonequilibrium systems at 
stationary state. This implies the following hypotheses: 

1. Tsallis entropy, as an information or uncertainty measure associated 
with the probability p^ is valid for some nonequilibrium nonextensive 
systems in steady or stationary states (time independent); 



2 



2. pi is the time invariant probability (or invariant measure[ll]) for the 
system under consideration to be found at the stationary state i; 

3. the Janyes' inference method (unbiased guess) of maximum entropy [12, 
13] applies for these states, i.e., the time independent probability distri- 
bution pi can be deduced from maximizing the entropy (uncertainty) 
under some constraints associated with, e.g., the normalization of Pi 
and the averages of physical quantities determining pj. In this work, p f 
is supposed to be uniquely determined by the energy Ei of the system 
at stationary state %. 

Due to the maximum Tsallis entropy, the framework of the conventional 
equilibrium thermodynamics can be copied for stationary nonequilibrium 
systems. Essential for this work is the establishment of zeroth law of ther- 
modynamics, i.e., the intensive variable (3 = ^f- (we call it inverse temper- 
ature from now on) is uniform everywhere at stationary state, where U q is 
the internal energy defined below in the section 3. It is worth noticing that, 
due to the different formalisms of NSM proposed in the past 15 years, the 
inverse temperature j3 has been given several definitions which sometimes 
cause confusion. A comment on this subject was given in ref. [14]. 

In what follows, we will focalize our attention on an important question 
about the validity of NSM for the systems containing subsystems with dif- 
ferent g's. Up to now, NSM has been applied only to composite systems 
containing subsystems with same q for which several versions of zeroth law 
of thermodynamics were established on the basis of Eq.(3)[15, 16, 17]. For 
the systems having different g's, the attempt to apply NSM [19, 20] using 
Eq.(3) is unfruitful since, firstly, no rigorous zeroth law can be formulated, 
secondly, no Tsallis entropy can be formulated as in Eq.(l) for a compound 
system A + B when A and B respectively have Tsallis entropy with their own 
qA and g#. This result implies that, in the actual formulation of NSM, there 
must be some fundamental aspects which are incomplete or handicapped. 
Nature contains nonextensive systems with different q values, as shown by 
the applications of NSM in the past years, and there must be also systems 
containing subsystems with different g's, such as a cluster of galaxies contain- 
ing different galaxies with different g's, a heterogeneous materials containing 
elements being physically and chemically different, a fractal containing two 
fractals of different dimensions and physics, or a social group including so- 
cially (culturally, economically or technically...) different subgroups, etc. If 
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NSM is useful only for the systems containing uniform distribution of q value, 
not only its application would be extremely limited, but worse, its general va- 
lidity would be fundamentally in danger[18] because, for different g-systems, 
no composability of physical quantity is possible, no measurement of tem- 
perature or pressure of a g-system is possible with ordinary thermometers or 
barometers which obey Boltzmann-Gibbs thermostatistics with q — 1 (or in 
order to measure the temperature of a g-system, you have to verify that the 
thermometer has the same q as the system!). This is equivalent to saying 
that NSM is useless in view of the diversity of physical or social systems. So 
this obstacle must be removed in order that NSM be a coherent and com- 
plete theory whose formalism, just like that of BGS, should remain the same 
whatever the level of the system, i.e., be hierarchically invariant. 

As a matter of fact, it is often forgotten that Eq.(3), a fundamental rule 
of NSM, was only established for the subsystems with same q. The afore- 
mentioned failure of the theory is in fact a failure of Eq.(3). The aim of this 
paper is to show that, by replacing Eq.(3) with a more general nonadditivity 
rule for Tsallis entropy, NSM can be applied to any system having different 
g's. We can in addition derive a calculus which uniquely determines the q of 
a composite systems from the g's of its subsystems. Without loss of general- 
ity, the discussion will be made with the g-exponential distribution given by 
Eq.(2). 

2 A more general pseudoadditivity for Tsallis 
entropy 

Now let us consider A and B, two nonextensive subsystems of a composite 
nonextensive system A + B. It has been proved[6] that the most general 
pseudoadditivity (or composability) of entropy or energy prescribed by zeroth 
law is the following : 

H[Q(A + B)] = H[Q{A)\ + H[Q(B)} + X Q H[Q(A)]H[Q(B)\, (5) 

which is in fact a very weak condition where H [Q] is just certain differentiable 
function satisfying H[0] = 0, Xq is a constant, and Q is either entropy S or 
internal energy U[7]. The functional form of H(Q) naturally depends on the 
nonadditive character of the system of interest. As shown in [7], for a given 
relationship Q = f(N) where N is the number of elements of the system, the 
finding of H(Q) is trivial. 
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As a matter of fact, Eq.(5) has been established [6, 7] for the class of sys- 
tems containing only subsystems having same q value. It must be generalized 
for the systems whose subsystems have different g's. This generalization is 
straightforward if we replace the Eq.(l) of reference [6], i.e., S(A + B) = 
f{S(A),S(B)} for uniform q, by H q [S q (A+B)] = f{H qA [S qA (A)], H qB [S qB (B)]} 

(or by H q [S q (A+B)] = H^[S u (A)]+H w [S^(B)]+g{H u [S u (A)lH m [S n (Bm 
where H q (S q ) is a functional depending on g's in the same way for the com- 
posite system as for the subsystems, where q, qA and q B are the parameters 
of the composite system A + B, the subsystems A and B, respectively. This 
functional H q (S q ) is necessary for the equality to hold in view of the different 
g's in the entropies of different subsystems. The function / (or g) is to be 
determined with the help of the zeroth law. Now repeating the mathematical 
treatments described in the references [6, 7], we find 

H q [Q(A + B)] = H qA [Q(A)\ + H qB [Q(B)] + X Q H qA [Q(A)]H qB [Q(B)}. (6) 

Eq. (5) turns out to be a special case of Eq. (6) for same q systems. In this case, 
if H q [S q ] is chosen to be identity function to give Eq.(3), then Tsallis entropy 
can be proved to be unique for this class of systems [4, 5] just as the unique- 
ness of Shannon information measure Si can be proved with the entropy 
additivity for independent subsystems having product joint probability [21]. 
This reasoning can be extended to the subsystems A and B each having its 
own q if one chooses H q (S q ) = J2T=iPl — 1 = (1 — q)S q . From Eq.(6), this 
choice implies : 

(l-q)S q (A + B) = (l-q A )S qA (A) + (l-q B )S qB (B) (7) 
+ X s (l-q A )(l-q B )S qA (A)S qB (B) 

where q, q A and q B are the parameters of the composite system A + B, the 
subsystems A and B, respectively. Eq.(7) is not derived from Eq.(5) but only 
postulated in order to generalize Eq. (3) according to the necessary condition 
Eq.(5) for the establishment of zeroth law. The uniqueness of Tsallis entropy 
for the systems obeying Eq.(7) can also be proved along the line of refer- 
ence [21] by replacing the Axiom 3 of the information additivity of [21] with 
Eq.(7) and the Eq.(19) of [5], i.e., [1 + (1 - q)L q (r m )] = [1 + (1 - q)L q (r)] m , 
by [1 + (1 - q c )L q (r m )} = UZi I 1 + i 1 ~ Qi) L M], where L q is the entropy 
functional, q c is the parameter of the composite system containing m sub- 
systems each having its own parameter qi (I = 1,2, ...m) and r(> 2) states. 
The details of this proof will be given in another paper. 



5 



Now let As = 1, Eq.(7) implies 



pt 3 (A + B)=p^(A) P r(B) (8) 

which can be called the generalized factorization of joint probability for the 
systems of different g's. A possible understanding of this relationship is that 
here the physical or effective probability is pj instead of Pi. This interpre- 
tation may shed light on why it has been necessary to use p\ or the escort 
probability['2] to calculated expectation. Eq.(8) can be written as the usual 
product probability p%j{A + B) — pi(A)pj(B) if and only if q = q A = qs for 
the systems having the same nonadditivity. 

Eq.(8) allows one to determine uniquely the parameter q for the composite 
system if q A , q B , Pi{A) and Pj(B) are given. By the normalization of the joint 
probability, we obtain the following relationship : 

W A W B 

Y,P" A/q (A)Y,pT /q (B) = l (9) 
i=l j=l 

which means q A < q < q B if q A < qs an d q = qA = 1b if Qa = 1b- In this 
way, for a composite system containing iV subsystems (k=l,2,...,N) having 
different qk, the parameter q is determined by 

N w k 

n e pZ'\a) = i, (io) 

k=l i k =l 

from which we can say that mini(qk) < q < maxi(qk)- 

Now if (A + B) is at stationary state, according to Jaynes principle of 
unbiased guess, dS(A + B) = 0. From Eq.(7), we get : 

(1 - q A )dS qA (A) | (1 - q B )dS QB (B) _ Q 



l + (l-q A )S gA (A) l + (l-q B )S qB (B) 

Eq.(ll) is essentially different from ^ + l+TT^fVe) = usec ^ fo r 

the subsystems with same q and will allow us to apply NSM to the systems 
in which the distribution of q value is not uniform. 



3 Zeroth law with nonadditive energy 

In general, for exact treatments of nonextensive systems within NSM, we 
have to deal with nonadditive energy that allows the existence of thermal 
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equilibrium and stationarity[7] and is compatible with the entropy nonad- 
ditivity. Usually, the nonadditivity of energy is determined by the product 
probability. Considering the distribution Eq. (2) and the generalized product 
probability Eq.(8), we get : 

hi[l-(l-q)P(A + B)E ij (A + B)} = (12) 



1-q 

Qa 



1-qA 



ln[l - (1 - q A )P(A)E i (A)} 



Qb 



\n[l-{l-q B )(3{B)E ] {B)]. 



Now let us define a deformed energy ef = nz^jjp m [l — (1 — Q)PEi] in order 
to write the distribution Eq.(2) in the usual exponential form 

Pi = ^exp[-(3e ( f ) ]. (13) 

The introduction of this distribution into the generalized product joint prob- 
ability Eq.(8) leads to : 

q/3(A + B)e[f(A + B) = q A f3(A)e ( f ) (A) + q B (3{B)ef {B). (14) 

Note that ef varies monotonically with so that dEi = leads to def = 
in accordance with the conservation of both E^ and ef\ 

In what follows, the zeroth law will be discussed within two possible for- 
malisms of NSM using the unnormalized expectation U q = J2i pfBi and the 
expectation given with escort probability U q = J2iPjBi/ Y^iPi [11]) respec- 
tively. In view of the generalized product joint probability Eq.(8) and the 
energy nonadditivity, the choice of these statistical expectations is necessary 
in order to split the composite expectations of the total system into those 
of the subsystems. The establishment of zeroth law would be impossible 
without this splitting. 



3.1 With unnormalized ^-expectation 

The formalism of NSM with unnormalized expectation U q = Y,i p\Bn was pro- 
posed by Tsallis and co-workers[l, 22] and has found many applications [23, 
24]. Recently, it has been shown [17] that, with this expectation, the inverse 
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temperature could be well defined by j3 = |^ for the systems of same q, so 
that all the details of the mathematics (including the statistical interpreta- 
tion of heat and work) of BGS remain the same within NSM. In what follows, 
it will be shown that the temperature can be defined between the subsystems 
with different q. 

The expectation of the deformed energy is given by u q = Y^iPi^l ■ m 
order to find a relationship between (3(A + B), (3(A) and (3(B), we assume 

qu q {A + B) qAU qA (A) q B u qB (B) 



EijPlM + B ) T^xp\ a {A) E J= iPf(B) 

which is consistent with Eq.(14) and the generalized product probability. It is 
in fact the analog of the additive energy in BGS. In general, this "additivity" 

does not means (3(A) = (3(B)] it only implies (3(A+B) = l ——^7 A n^/~mf ^ 

where u' g = qAU qA / Yl,i=\P\ A ■ From energy conservation, we obtain du q (A + 
B) = and 

q A du qA (A) q B du qB (B) 

ztipr(A) + E^pf(s) • 1 ) 

Now in order to determine (3 in the distribution of Eq.(2) or Eq.(13), we 
take a deformed entropy s q defined by Taneja[25, 26] : 

10 

S q = -J^Pi lli Pi ( 17 ) 

1=1 

but here Pi is given by Eq.(13). This deformed entropy is concave and has 
the following nonadditivity : 

qs q (A + B) qAS qA (A) q B s qB (B) 

Z%pU a + B ) ^pV(A) + EiPfdJ) 1 } 

which leads to, at maximum entropy, I^~%a^ + = 0- Combining 

this with Eq.(16) and considering s q = J2iPi hi Z q + (3u q , we get 

(3(A) = P(B) (19) 

where the inverse temperature is given by (3 = ^ = It is worth 

noticing that the equality of the temperatures of A and B at stationary state 
is independent of whether q A and q B are equal. 



S 



3.2 With escort probability 

The above description of thermal equilibrium can also be made with escort 
probability proposed in [11]. Let us define a deformed average energy with 
escort probability fi q = J2iPi e i^ I JliPi, where pi is the probability given by 
Eq.(2). Still assuming qfi q (A + B) = qAfi qA (A) + qBH qB (B) which means, 
from the nonadditivity Eq.(14), (3(A + B) = g^gg^g^H , one 
gets 

q A dfi qA (A) + q B dfi qB {B) = 0. (20) 

in accordance with the conservation of ef or fi q . 

Now let us take the entropy ip q defined by Aczel and Daroczy[25, 27] as 
a deformation of the nonextensive entropy S q : 

1 w 

^ = -v5-3E^ft (21) 
2^i=lPi i=\ 

which naturally arises from the escort probability. ip q can be maximized with 
appropriate constraints associated with the deformed average energy fi q to 
give the distribution of Eq.(13). 

The above deformed entropy has the following nonadditivity : 

<p/> q (A + B) = q A if> qA (A) + QB^ qB (B), (22) 

so that at maximum entropy dip q (A+B) = 0, we have qAdip qA (A)+q B dip qB (B) ■ 
0. Combining this with equation Eq.(20), one gets 

d^jA) = djj QB (B) 

d^ A (A) dfi qB (BY 1 ) 

On the other hand, from Eqs.(13) and (21), it is easy to write 

i; q = \nZ q + P[i q (24) 
which means j3 = ^f-. We finally have 

/3(A) = (3(B) (25) 
which is also independent of whether q A and q B are equal. 
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4 An alternative method 



In the above discussions, the deformed entropies s q and ip q are introduced in 
order to discuss nonextensive statistics with the language of extensive statis- 
tics BGS. s q and ip q are in fact alternative expressions of S q associated with 
the deformed energy u q and \i q , respectively. They can be maximized, just 
as S q , to give g-exponential distribution [28] in terms of the deformed energy 
e^f 1 [Eq.(13)]. The description of the stationarity characterized by Eq.(19) 
or Eq.(25) is of course independent of their introduction. In what follows, I 
we briefly present an alternative method without these deformations. 

Using the product joint probability Eq.(8) and the relationship YUpI = 
Z\~ q + (1 — q)ftU q calculated from the distribution Eq.(2) and the unnormal- 
ized expectation U q , we get 

Z\~\A + B) + (1 - q)(3(A + B)U q (A + B) (26) 
= [Z]^{A) + (1 - q A )(3(A)U qA (A)][Zl^(B) + (1 - q B )/3(B)U qg (B)]. 

Considering the total energy conservation dU q (A + B) = 0, we obtain 

(l-q A )(3(A)dU qA (A) , (l-q B )P(B)dU qB (B) 



which suggests following energy nonadditivity 

(1 - q A )dU qA (A) (1 - q B )dU qB (B) 
H,pi A {A) + Z iP t B (B) 



(27) 



(28) 



as the analog of the additive energy dU(A)+dU (B) = of BGS. From Eq.(28) 
and Eq.(ll) follows (3(A) = (3{B) with (3(A) = and (3(B) = g^gg . 



5 Conclusion 

We have discussed the establishment of the zeroth law for stationary state 
between nonextensive systems with different g's. According to our starting 
hypotheses, this stationary state maximizes Tsallis entropy so that the zeroth 
law of the equilibrium thermodynamics can apply. This work is carried out 
with a generalized nonadditivity rule of Tsallis entropy which reduces to the 
usual one if all the subsystems have the same q. It should be noticed that 
the generalized product joint probability Eq.(8) derived from the generalized 
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nonadditivity of entropy suggests two possible expectation calculus of NSM, 
the unnormalized g-expectation and the escort probability g-expectation, in 
order to establish zeroth law. A main result of this work is to show that 
the intensive variable (3 is uniform in the nonextensive systems in stationary 
state, independently of whether or not the systems contain subsystems with 
different q values. So NSM theory based on the maximization of Tsallis 
entropy is, contrary to some belief, also valid for the systems containing non 
uniform distribution of q values. 
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